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1^-^ , Abstract. We present the results of N2H+(l-0) observations of 35 dense molecular cloud cores from the northern and southern 

' hemispheres where massive stars and star clusters are formed. Line emission has been detected in 33 sources, for 28 sources 

I detailed maps have been obtained. Peak N2H+ column densities lie in the range: 3.6 • 10'' - 1.5 • lO''' cm"^. Intensity ratios 

7—1 . of (01-12) to (23-12) hyperfine components are slightly higher than the LTE value. The optical depth of (23-12) component 

^ ' toward peak intensity positions of 10 sources is ~ 0.2 - 1. In many cases the cores have elongated or more complex structures 

0^ , with several emission peaks. In total, 47 clumps have been revealed in 26 sources. Their sizes lie in the range 0.3-2.1 pc, 

the range of virial masses is ~ 30 — 3000 Mq. Mean N2H^ abundance for 36 clumps is 5 ■ 10 Integrated intensity maps 
with aspect ratios < 2 have been fitted with a power-law radial distribution ~ r^'' convolved with the telescope beam. Mean 
power-law index for 25 clumps is close to 1.3. For reduced maps where positions of low intensity are rejected mean power-law 
index is close to unity corresponding to the ~ r^^ density profile provided N2H+ excitation conditions do not vary inside these 
I regions. In those cases where we have relatively extensive and high quality maps, line widths of the cores either decrease or stay 

— , ■ constant with distance from the center, implying an enhanced dynamical activity in the center There is a correlation between 

i-C ' total velocity gradient direction and elongation angle of the cores. However, the ratio of rotational to gravitational energy is 

too low (4-10 ''-7. MO ^) for rotation to play a significant role in the dynamics of the cores. A correlation between mean line 
^ ■ widths and sizes of clumps has been found. A comparison with physical parameters of low-mass cores is given. 

^: 
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1. Introduction anomalies in HCN (Pirogov [1999) and ammonia (Stutzki & 
^ Winnewisser 1985 ) in warm clouds. The degree of small-scale 
. It is well-known that stars in some cases are formed in clus- clumpiness is probably related to the number of newborn stars. 
- - ' ters and isolation in others, however, the reasons of such a di- 
chotomy are far from being understood at present. The number Important parameters for describing evolution of star form- 
of stars depends apparently on physical parameters of dense ing cores in modern theoretical models are radial dependences 
molecular cloud cores where stars are formed. It is known, in of density and velocity dispersion in envelopes that surround 
particular, that clusters are formed in the most massive and tur- newborn stars. A popular model for massive star formation is 
bulent cores. High-mass stars are encountered solely in such based on the logatropic equation of state which incorporates 
cores. An important feature of massive cores is their clumpy non-thermal gas motions (McLaughhn & Pudritz[T996j. It re- 
structure which extends from large scales down to the small- suits in the ~ density distribution. Velocity dispersion in 
est ones even unresolved by present-day telescopes. An exis- this model is nearly constant in the center and rises with dis- 
tence of small-scale clumpiness is implied, for example, by 1- tance in the outer regions as ~ r^'\ In contrast, the standard 
2 orders of magnitude difference between mean densities and theory of low-mass star formation from a singular isothermal 
densi ties de rived from multit ransitional analys is (e g . Bergin sphere (Shu 1977 ) predicts the ~ density profile in the static 
et al. 119961 Zinchenko 119951 Zinchenko et al. I1998II: by the outer envelope. Both models predict ~ r^'^ and ~ r '^^ density 
absence of significant broadening of presumably very opti- and velocity profiles, respectively, for the regions involved in 
cally thick CS lines as compared to the rarer C^'^S isotopomer free-falling collapse (Shu .1977. McLaughhn & Pudi-itz [T997) . 
lines (Zinchenko et al. 1998 ); by the hyperfine intensity ratios Therefore, studies of density and velocity dispersion profiles in 

the cores are essential for discrimination between different the- 

Send ojfprint requests to: L. Pirogov oretical models and for determination of their dynamical status. 
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Recently, density structure of large samples of objects have 
been explored in several studies. Analysis of dust continuum 
emission in low-mass cores where stars are formed mostly in 
isolation, reveals a density power-law index of or ~ 1 .5 - 2 (e.g. 
Motte & Andre IMm Shirley et al. f30D3| being more or less 
in accordance with the isothermal sphere model. Dust contin- 
uum emission studies in high-mass star forming regions have 
resulted in the following power-law indices: {a} = 1.6(0.5) 
(Beuther et al. 12003 and (a) = 1.8(0.4) (Mueller et al. 2002 1, 
while molecular multi-line modeling of several massive star 
forming regions (Van der Tak et al. l2000> gave somewhat lower 
values: a = 1.0-1.5. 

In their study of velocity dispersion profiles in low-mass 
cores, Goodman et al. ( I1998> found nearly constant velocity 
dispersion within so-called "coherent regions" with sizes of 
~ 0.1 pc and a rise of velocity dispersion in the outer re- 
gions with power-law index ^ 0.2, which is more consistent 
with the logatropic model. In contrast, molecular line studies of 
high-mass cores show highly non-thermal lines of decreasing 
widths with distance from center (Zinchenko l 1 9951 Lapinov et 
al. 119981 Fontani et al. 2002 1, which apart from optical depth 
effects could imply a higher degree of gas dynamical activity 
in central regions. 

Recently, about 60 nearby low-mass cores were mapped in 
the N2H+(l-0) line (CaselU et al. 2002 hereafter Paper I). As a 
result, a number of physical parameters including density and 
velocity dispersion profiles were derived. N2H+ was chosen as 
a tracer of quiescent high density gas (critical density for exci- 
tation of the 7=1-0 transition at a kinetic temperature of 40 K 
is 4.4 ■ 10^ cm"^). An advantage of the N2H+(l-0) transition is 
the hyperfine structure which enables a reliable determination 
of the optical depth and improves estimates of related physical 
parameters of the dense gas that surrounds young stars. It is 
known that the N2H^(l-0) lines have relatively strong intensi- 
ties and moderate (t ~ 1) optical depths both in cold and warm 
clouds (Womack et al. "1992) associated with low-mass and 
high-mass star forming regions, respectively. The goal of the 
present study is to derive physical parameters for a represen- 
tative sample of dense cores in high-mass star forming regions 
from observations in the N2H^(l-0) line. The derived param- 
eters will be compared with those of low-mass cores (Paper I) 
to help to understand the details of the star formation process. 

A sample list for the study has been taken from the ammo- 
nia database compiled by P. Myers and co-workers (Jijina et 
al. 1999) where data on cluster existence is included. This list 
overlaps with a sample of massive cores earlier mapped in CS, 
both in the northern and southern hemispheres (Zinchenko et 
al. [T994|[T995 'T998' JuvelalW). In the present paper we give 
the results of N2H^(l-0) observations towards a sample of 35 
massive dense cores where star formation in clusters is going 
on. Their physical parameters including density and velocity 
dispersion profiles are derived and compared with parameters 
of low-mass cores (Paper I). 

2. Observations 

The N2H^(l-0) observations at the frequency 93.2 GHz have 
been done in April, 2000 at 20-m OSO radiotelescope and in 



February, 2001 at 15-m SEST radiotelescope. The observing 
parameters are given in Table[2 

For OSO we took into account beam efficiency dependence 
on source elevation and give an average value. We also give 
the ranges of system temperatures which depend on source el- 
evation and weather conditions. At OSO we used both a fil- 
terbank spectrum analyzer (256 channels) and autocorrelator 
(1600 channels) while at SEST an acousto-optical spectrum 
analyzer (1000 channels) was used. Pointing was regularly 
checked by SiO maser observations and usually was better than 
5". 

Table 1. N2H^(l-0) observing parameters. 



Telescope 


A0 


'7mb 




SV 




(") 




(K) 


(km s-') 


OSO-20m 


40 


~ 0.5 


250-500 


0.8; 0.16 


SEST- 15m 


55 


0.75 


200-300 


0.14 



The objects for observations were selected from the list of 
dense cores observed earlier in the CS(2-1) line (Zinchenko et 
al. lT994|[T995l[T998l Juvela [T996l as well as from the ammonia 
database (Jijina et al. 1 19991 1 according to the following criteria: 
an existence of embedded clusters of stars detected in infrared 
(Jijina et al. il999i Hodapp IT^ Carpenter et al. il993< . dis- 
tances to the objects not far than 5 kpc and an existence of 
extended CS or NH3 emission regions (> 2') in order to get de- 
tailed N2H^(l-0) maps. In the cases of no data on the clusters 
existence, the sources with high IR-luminosities were selected 
(L > 10^ Lq), which can be an indirect indication of cluster 
(Jijina et al. ll999j . In total, 35 objects that satisfy these criteria 
were selected. The source list with coordinates and distances is 
given in Table|2] The number of points per map and the names 
of associated objects are also given. Mapping has been done 
with 20" grid spacing (Nyquist sampled maps) with the excep- 
tion of S 187, G 133.69 and S 153 for which the spacing was 
40" (beam sampled maps). 

3. Results 

3.1. Integrated intensity maps 

The N2H^(l-0) emission has been detected in 33 of 35 sources 
(except LkHalOl and LkHQ'25). Maps of intensities integrated 
over all hyperfine components for 31 object are given in Fig.Q] 
The maps of AFGL 490 and G 142.00 (which actually be- 
long to a single molecular complex), S 88B, G 267.94 and 
G 305.36 have not been completed (the latter two are not shown 
in Fig.Q. The maps of S 187, G 133.69 and AFGL 6366 each 
consist of two separated clumps at the 30% peak intensity level. 
Some others (e.g. S 255, G 316.77, G 345.41, G 351.41) con- 
tain intensity peaks separated at the half maximum intensity 
level. The maps of RNO IB, NGC 281, G 173.48, S 76E, 
G 270.26 indicate spherical symmetry while others demon- 
strate elongated, curved or more complex structures. Per 4 and 
L 1251T4 show the most narrow and weak lines and actu- 
ally belong to low-mass star forming cores (see Paper I). The 
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Table 2. Source list. 



Source 


KA (IVjUJ 


uec {lyjv) 


u 


Points 


Associated 




(( ) ( ) ( ) ) 


( ) 


(kpc) 


per map 


objects 


RNO IB 


00 33 53.3 


63 12 32 


0.85" 


70 


L1287, G 121.30+0.66, IRAS 00338+6312 


NGC 281 


00 49 29.2 


56 17 36 


3.5* 


34 


S184, G 123.07-6.31 


S 187 


01 19 58.0 


61 34 32 


1.0 


67 


G 126.68-0.83 


G 133.69+1.22 


02 21 40.8 


61 53 26 


2.1 


52 


W3 


S 199 


02 57 36.9 


60 17 32 


2.1 


46 


IC1848, AFGL 4029 


S201 


02 59 22.4 


60 16 12 


2.1 ' 


41 


AFGL 437, G 138.50+1.64 


AFGL 490 


03 23 34.0 


58 33 47 


0.9 f 


18 




G142.00+1.83 


03 23 41.0 


58 36 52 


0.9 f 


37 


AFGL 490 


Per 4 


03 26 12.1 


31 17 13 


0.35 s 


78 




LkHalOl 


04 26 57.3 


35 09 56 


0.8'' 


9 


AFGL 585, S222 


G 173.48+2.45 


05 35 51.3 


35 44 16 


2.3'' 


54 


S231 


G 173.58+2.44 


05 36 06.3 


35 39 21 


2.3'' 


41 




AFGL 6366 


06 05 40.1 


21 31 16 


2.0' 


50 


S247 


S255 


06 09 58.2 


18 00 17 


2.5'' 


67 


AFGL 896 


LkHQ'25 


06 37 59.5 


09 50 53 


0.8^ 


1 


AFGL 4519S, NGC 2264, S273 


S76E 


18 53 47.0 


07 49 25 


2.1*= 


61 


G40.50+2.54 


S 88B 


19 44 41.4 


25 05 17 


2.0'' 


25 




DR21 


20 37 13.5 


42 12 15 


3.0' 


56 


W75 


S 140 


22 17 41.1 


63 03 44 


0.9'' 


82 


LI 204 


L 125 1T4 


22 27 10.2 


74 58 00 


0.35"' 


119 




S 153 


22 56 42.0 


58 28 45 


4.0'' 


42 




G 264.28+1.48 


08 54 39.0 


-42 53 30 


1.4" 


43 




G 265.14+1.45 


08 57 36.3 


-43 33 38 


1.7" 


64 




G 267.94-1.06 


08 57 21.7 


-47 19 04 


0.7" 


9 




G 268.42-0.85 


09 00 12.1 


-47 32 07 


1.3" 


60 




G 269.11-1.12 


09 01 51.6 


-48 16 42 


2.6" 


53 




G 270.26+0.83 


09 14 58.0 


-47 44 00 


2.6" 


27 




G 285.26-0.05 


10 29 36.8 


-57 46 40 


4.7" 


59 




G 291.27-0.71 


1 1 09 42.0 


-61 01 55 


2.7° 


104 




G 294.97-1.73 


11 36 51.6 


-63 12 09 


1.2" 


39 




G 305.36+0.15 


13 09 21.0 


-62 21 43 


4.2" 


18 




G 316.77-0.02 


14 41 10.4 


-59 35 50 


3.1" 


136 




G 345.01+1.80 


16 53 17.3 


-40 09 23 


2.1" 


103 


RCW 116 


G 345.41-0.94 


17 06 02.2 


-41 32 06 


2.8" 


71 


RCW 117 


G 351.41+0.64 


17 17 32.3 


-35 44 03 


1.7" 


88 


NGC 63341 



" Yang et al. fTMD , ' Martin & Henning HM^ . " Fich & Blitz JM?), Blitz et al. assumed, f Harvey « Paper I, Lada 

(T98?l , ' Carpenter et al. IT9931 , ' Snell et al. (T990l . * Plume et al. (1992}, ' Harvey et al. Il977l . Toth & Walmsley (1995), " Zinchenko et 

al. (1995}, " Brand & Blitz (1993), " Juvela (T595), 1 Neckel jT978j 



clumpy structures visible in the maps of these two sources (see 
Fig.[0 are primarily caused by noise. The OSO and SEST beam 
sizes are shown in Fig.^on RNO IB and G 268.42 maps, re- 
spectively. 

Most of the sources have associated IRAS point sources. 
The positions of IRAS sources are indicated by stars and the 
uncertainty ellipses corresponding to 95% confidence level are 
also shown. In 17 objects IRAS sources lie within the half 
maximum level (marked by thick black lines) or coincide with 
a secondary intensity peak. Although there is no IRAS point 
source in DR 21, it is a well-known star forming molecu- 
lar cloud with compact H II regions and strong extended far- 
infrared emission (Harvey et al. I1979> which is close in shape 
to the N2H^(l-0) integrated intensity distribution. 

Peak integrated intensities are given in Table |3 For the 
sources having separate clumps at the 30% level, peak intensity 



for each clump is given. The uncertainties of integrated inten- 
sities have been calculated as AFmb VA^^VcH, where ArMB 
- noise level in channel without line (after baseline subtract- 
ing), A^ch - number of line channels, 5Vch - velocity resolution. 
For the sources observed at OSO we used high-resolution au- 
tocorrelator data, with the exception of S 88B for which only 
filterbank data are available. 

3.2. N2H'^(1-0) line parameters and peak column 
densities 

N2H^(l-0) spectra consist of seven hyperfine components 
(Womack et al. 119921 Caselli et al. 1995. 1. However, due 
to overlap of closely spaced components only three lines 
(two of which are actually triplets) have been observed in 
most cases (with exception of LI 251 and Per4). As an ex- 
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Table 3. N2H^(l-0) line parameters and LTE column densities 



Source (Aff", A(5") / rMB(23-12) Rn R02 Klsr(23-12) AV yV(N2H+) 

(Kkms-') (K) (kms-') (km s"') (cm-^) 

(10'2) 



RNO IB (0,0) 


27.0(0.4) 


2.88(0.04) 


0.66(0.01) 


0.47(0.02) 


-17.55(0.01) 


2.39(0.03) 


31.1 


NGC 281 (0,0) 


27.2(0.4) 


2.29(0.03) 


0.68(0.02) 


0.48(0.02) 


-30.36(0.02) 


2.91(0.04) 


31.2 


S 187 (0,0) 


6.3(0.4) 


1.6(0.1) 


0.68(0.07) 


0.57(0.08) 


-13.33(0.03) 


0.91(0.06) 


7.2 


S 187 (120,-80) 


4.9(0.4) 


0.73(0.07) 


0.7(0.1) 


0.6(0.1) 


-13.59(0.05) 


1.5(0.2) 


5.6 


G 133.69+1.22 (160,-160) 


11.7(0.7) 


1.4(0.1) 


0.62(0.07) 


0.6(0.1) 


-38.38(0.05) 


1.8(0.1) 


13.4 


G 133.69+1.22 (40,-80) 


8.0(0.7) 


0.38(0.05) 


1.0(0.2) 


0.8(0.2) 


-41.9(0.2) 


3.6(0.4) 


9.2 


S 199 (0,-20) 


9.7(0.6) 


0.94(0.06) 


0.64(0.07) 


0.6(0.1) 


-38.02(0.06) 


2.3(0.2) 


11.2 


S 201 (40,40) 


3.9(0.4) 


0.26(0.03) 


1.2(0.2) 


0.8(0.2) 


-37.3(0.1) 


2.8(0.3) 


4.5 


AFGL 490 (-20,0) 


7.6(0.6) 


1.5(0.1) 


0.67(0.09) 


0.6(0.1) 


-12.24(0.04) 


1.1(0.1) 


8.7 


G142.00 (-20,-80) 


10.9(0.7) 


1.2(0.1) 


0.7(0.1) 


0.4(0.1) 


-13.00(0.06) 


1.8(0.2) 


12.5 


Per 4 (0,0) 


4.5(0.3) 


2.07(0.09) 


0.69(0.05) 


0.58(0.05) 


7.62(0.01) 


0.48(0.02) 


5.2 


G 173.48+2.45 (0,0) 


43.4(0.5) 


3.78(0.04) 


0.72(0.01) 


0.53(0.02) 


-16.12(0.01) 


2.70(0.03) 


49.9 


G 173.58+2.44 (0,-20) 


12.2(0.4) 


1.71(0.06) 


0.74(0.04) 


0.56(0.05) 


-16.42(0.02) 


1.63(0.07) 


14.0 


AFGL 6366 (0,0) 


14.7(0.3) 


1.56(0.03) 


0.70(0.02) 


0.56(0.03) 


2.68(0.02) 


2.16(0.05) 


16.9 


AFGL 6366 (80,60) 


10.5(0.4) 


1.45(0.06) 


0.67(0.05) 


0.53(0.06) 


2.10(0.03) 


1.63(0.08) 


12.0 


S 255 (0,60) 


27.5(0.6) 


2.58(0.05) 


0.67(0.02) 


0.53(0.03) 


8.97(0.02) 


2.61(0.06) 


31.6 


S 76E (-20,-20) 


79.3(0.7) 


6.03(0.05) 


0.70(0.01) 


0.53(0.02) 


32.44(0.01) 


3.10(0.03) 


91.1 


S 88B (20,20) 


4.3(0.8) 


0.5(0.2) 


0.6(0.4) 


0.7(0.4) 


22.7(0.1) 


1.3(0.4) 


5.6 


DR 21 (0,^0) 


106.3(0.7) 


6.93(0.04) 


0.72(0.01) 


0.66(0.01) 


-3.25(0.01) 


3.51(0.02) 


122.0 


S 140 (60,20) 


48.2(0.7) 


6.53(0.09) 


0.64(0.02) 


0.56(0.02) 


-7.03(0.01) 


1.73(0.03) 


55.4 


L 125 1T4 (-60,20) 


3.1(0.3) 


1.0(0.1) 


0.65(0.13) 


0.37(0.13) 


-3.81(0.03) 


0.49(0.05) 


3.6 


S 153 (0,40) 


21.2(1.0) 


2.0(0.1) 


0.76(0.06) 


0.55(0.08) 


-51.23(0.05) 


2.35(0.12) 


24.4 


G 264.28+1.48 (20,-20) 


6.2(0.1) 


0.64(0.01) 


0.70(0.02) 


0.51(0.02) 


6.64(0.01) 


2.17(0.04) 


7.1 


G 265.14+1.45 (0,-60) 


15.8(0.2) 


1.61(0.02) 


0.75(0.01) 


0.57(0.02) 


7.40(0.01) 


2.11(0.03) 


18.2 


G 267.94-1.06 (0,20) 


3.4(0.2) 


0.56(0.04) 


0.53(0.07) 


0.37(0.08) 


3.29(0.04) 


1.36(0.12) 


3.9 


G 268.42-0.85 (0,-20) 


10.9(0.2) 


1.13(0.02) 


0.66(0.02) 


0.44(0.03) 


3.40(0.02) 


2.56(0.06) 


12.5 


G 269.11-1.12 (0,40) 


24.3(0.2) 


1.66(0.01) 


0.76(0.01) 


0.56(0.02) 


10.49(0.01) 


3.23(0.03) 


28.0 


G 270.26+0.83 (-20,20) 


20.9(0.2) 


1.44(0.02) 


0.73(0.01) 


0.49(0.02) 


9.69(0.02) 


3.42(0.03) 


24.0 


G 285.26-0.05 (-60,-60) 


12.4(0.2) 


0.92(0.01) 


0.63(0.02) 


0.55(0.02) 


2.54(0.02) 


3.12(0.04) 


14.1 


G 291.27-0.71 (-80,-100) 


36.1(0.2) 


2.91(0.01) 


0.69(0.01) 


0.53(0.01) 


-24.36(0.01) 


3.00(0.02) 


41.6 


G 294.97-1.73 (20,-20) 


15.6(0.2) 


1.53(0.02) 


0.69(0.01) 


0.41(0.02) 


-8.05(0.01) 


2.53(0.03) 


17.9 


G 305.36+0.15 (20,40) 


7.0(0.2) 


0.47(0.02) 


0.69(0.05) 


0.64(0.07) 


-38.54(0.05) 


3.2(0.1) 


8.1 


G 316.77-0.02 (0,20) 


51.6(0.2) 


3.53(0.01) 


0.72(0.01) 


0.60(0.01) 


-39.56(0.01) 


3.39(0.01) 


59.3 


G 345.01 + 1.80 (-20,20) 


58.5(0.2) 


4.07(0.02) 


0.65(0.01) 


0.56(0.01) 


-14.08(0.01) 


3.55(0.01) 


67.3 


G 345.41-0.94 (80,-20) 


13.7(0.2) 


1.05(0.01) 


0.64(0.01) 


0.47(0.02) 


-22.31(0.02) 


3.25(0.04) 


15.7 


G 351.41+0.64 (20,100) 


128.7(0.2) 


6.35(0.03) 


0.78(0.01) 


0.79(0.01) 


^.23(0.01) 


4.38(0.02) 


148.0 



ample, the N2H^(l-0) spectrum toward the peak position in 
G 351.41+0.64 is shown in Fig. |2] with hyperfine compo- 
nents indicated. The relative amplitudes of the components 
correspond to the LTE optically thin case. Each component is 
marked by two pairs of numbers {F[ F' - Fi F) corresponding 
to initial and final states of the transition (Caselli et al. ll995> . 

After removing baselines from the data (low-order poly- 
nomials) we performed gaussian fitting. While choosing three 
gaussians for fitting the data can give only crude informa- 
tion about position and amplitudes of the components, fitting 
by seven independent gaussians with fixed separations usually 
leads to incorrect results (negative amplitudes for some of the 
overlapped components). We decided to fit the data by seven 
gaussians, only three of which having independent amplitudes, 
whereas the amplitudes of the overlapping components have 
been assumed proportional to their statistical weights. This 



method proved to be useful for treating both LTE and non-LTE 
cases and gave lowest values in most cases. 

The derived parameters are (Table O: main beam bright- 
ness temperature of the (23-12) component, peak intensity ra- 
tios of the (12-12) to (23-12) and (01-12) to (23-12) transi- 
tions {R\2 and 7?()2, respectively), velocity of the (23-12) com- 
ponent and line width at half maximum level, considered equal 
for all components. For most of the sources the R\2 ratios 
do not significantly differ from the optically thin value (5:7) 
while the i?02 ratios are usually higher than expected in the 
optically thin case (3:7). Observed R12 versus 7?o2 ratios to- 
gether with the LTE curve are shown in Fig.|3 The mean ratios 
for the sources given in Table|3]are: {Rn) - 0.70(0.01) and 
(^02) = 0.57(0.03). Similar ratios apply for each source when 
all map positions with significant emission are averaged, i.e., 
{Rn) is close to 5:7 while {Rqi) is higher than 3:7. Fitting by 
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Fig. 1. N2H^(l-0) integrated intensity maps. Intensity contour 
ranges from 30% to 90% of peak intensity, thick contours 
correspond to half maximum levels. Observed positions are 
marked by dots and IRAS point sources are marked by stars. 
The uncertainty ellipses corresponding to 95% confidence level 
in IRAS point source position are also shown 

three independent gaussians also reveals an enhanced intensity 
of the (01-12) component compared with the LTE case. 

No evidence of broad line wings as indicators of outflow 
activity has been found in the spectra. The N2H^(l-0) spectra 
towards several positions of G 316.77, G 345.01 and G 351.41 
are highly non-gaussian apparently associated with more than 
one clump along the line of sight. Fitting these spectra by the 
procedure described above resulted in much higher line widths 
than for the rest of spectra. No attempts have been made to 
decompose these spectra into several velocity components due 
to the ambiguity of such a procedure, and we have excluded 
these sources from our velocity analysis. 

In order to derive line optical depths in the case of closely 
spaced and overlapping components the following function 
could be used: 

T(v) = [^7(rEx,)^^^-/(rBc)]-[l-exp(-^T,(v))],(l) 

where J(T) = ^[exp(|^) - l]"', Tbg = 2.7 K is the cosmic 
background radiation temperature, Tex; and t,(v) are excita- 




AaC) 



Fig. 1. continued 



tion temperature and optical depth of the i-th component at the 
frequency v, respectively, and summing is going over all com- 
ponents. For the case of the observed N2H^(l-0) spectra the 
fitting function Q could be simplified. It is known that exci- 
tation temperatures of the overlapping lines always tend to be 
equal due to photon exchanges (e.g. Guilloteau et al. I198I> . 
Assuming equal excitation temperatures for the overlapping 
components within triplets, and, that the (01-12) line and the 
triplets do not overlap, the fitted function becomes: 



T(y) = J(Tex,) - J(Tbc)] ■ [l - exp(- ^ t^v))] 



(2) 



where the first sum is going over non-overlapping groups of 
components and the second is going over overlapping compo- 
nents within individual group. Tit{v) is the optical depth of the 
i-th group and k-th overlapping component at the frequency v. 
Thus, the fitted parameters are: three excitation temperatures, 
optical depth of the central (23-12) component as well as line 
width and velocity of the central component. 

The spectra toward the positions given in Table |3] as well 
as toward several additional emission peaks have been fitted by 
the function Q. In three sources with low S/N ratios (G 133.69, 
S 88B and L 1251) the solutions did not converge. In 10 
cases the derived optical depth exceeds 3cr. For the remain- 
ing sources an upper 3cr level has been derived and it exceeds 
unity in 8 cases (in the case of S 201 t(23-12)<3, the high- 
est upper limit among the sources). The results of fitting for 
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10 sources (excitation temperatures, velocities, line widths and 
optical depths) are given in Table|4] The optical depth and exci- 
tation temperature ranges for different sources are: ~ 0.2- 1 and 
~ 7 - 29 K, respectively. The difference in excitation temper- 
ature between different groups of components is < 5.6 K and 
in several cases 7'ex(01-12) is the highest one. As the sources 
given in Table |3 are among those with highest intensities we 
conclude that Tex ^ 25 K and t(23 - 12) ^ 1 for the studied 
sample. It should be emphasized that we have assumed a beam 
filling factor of 1 in all cases; if smaller, excitation temperatures 
are higher. 

Since non-LTE effects seem small, we have assumed LTE 
and the optically thin approximation to calculate N2H+ col- 
umn densities, A^(N2H^), for all sources in Table |3] Excitation 
temperatures were taken equal to 10 K for all sources, the 
value at which LTE column densities are close to minimum. 
At Tex = 20 K column densities are 1 .6 times higher The cal- 
culated values, which could be treated as lower limits, are given 
in the last column of Table|3] 

3.3. Clump deconvolution 

Most of the maps in Fig. [2 show elongated or more complex 
structures with prominent secondary peaks. In order to esti- 
mate core sizes, two-dimensional (2D) gaussian fittings have 
been used. As a first step we tried to fit simultaneously a num- 
ber of circular 2D gaussians to each map in order to reveal dis- 



tinct clumps. The initial guesses of number of components and 
their locations were simply based on inspections of the maps 
in Fig. [2 and strip scans, obtained where required. In the sec- 
ond step the circular gaussians were replaced by elliptical ones, 
adding two more parameters to be determined: axial ratio and 
position angle. Finally, core sizes were estimated by decon- 
volving the geometric mean of the extents of the fitted ellip- 
tical gaussians with the beam size, i.e., C = VG2 - B2 where 
C, G, and B are the full half widths of the core, observed geo- 
metric mean of the extents and beam, respectively. In general, 
1-2 clumps per map were found but in several sources (S 187, 
G 291.27, G 316.77, G 345.01 and G 351.41) three clumps 
have been revealed. We have excluded Per 4 and L 1251T4 
from this analysis, (the former has already been analyzed in 
Paper I) as well as sources with small (< 25 points per map) 
and incomplete maps. In total, 47 clumps in 26 objects have 
been revealed. The relative coordinates of the clumps, their ax- 
ial ratios, angular and linear sizes calculated at half maximum 
level are given in Table |5l together with fitting uncertainties. 
We have not taken into account possible uncertainties in source 
distances while calculating linear sizes. Clumps in the same 
core are marked by numbers. Angular sizes of some clumps 
are smaller than the telescope beam. These clumps and those 
having high uncertainties in size have been excluded from cal- 
culations of physical parameters (Section [3.4> . For four cores 
where all revealed clumps are smaller than the telescope beam 
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Table 4. N2H^(l-0) line parameters with optical depths 



Source (Aor", A5") 


rEx(12-12) 


rEx(23-12) 


rEx(01-12) 


Vlsr(23-12) 


AV 


r(23-12) 




(K) 


(K) 


(K) 


(km s"') 


(km s ') 




RNO IB (0,0) 


7.75(0.37) 


9.18(0.38) 


7.32(0.50) 


-17.54(0.01) 


1.94(0.06) 


1.0(0.1) 


S 255 (0,60) 


9.3(1.6) 


10.6(1.6) 


10.1(2.2) 


8.96(0.02) 


2.31(0.10) 


0.5(0.2) 


S 76E (-20,-20) 


24.4(4.2) 


27.0(4.3) 


27.1(5.4) 


32.44(0.01) 


2.86(0.06) 


0.3(0.1) 


DR21 (0,^0) 


20.3(1.3) 


22.8(1.3) 


25.0(2.1) 


-3.25(0.01) 


3.09(0.04) 


0.6(0.1) 


S 140 (60,20) 


21.4(3.7) 


25.4(4.1) 


27.0(5.5) 


-7.03(0.01) 


1.62(0.05) 


0.4(0.1) 


G 270.26+0.83 (-20,20) 


8.8(1.5) 


9.2(1.5) 


8.8(1.8) 


9.69(0.02) 


3.17(0.08) 


0.3(0.1) 


G 291.27-0.71 (-80,-100) 


24.2(6.0) 


25.9(6.2) 


29.1(7.7) 


-24.36(0.01) 


2.90(0.03) 


0.15(0.05) 


G 316.77-0.02 (240,-60) 


8.9(0.6) 


10.1(0.6) 


10.1(0.9) 


-37.94(0.01) 


3.61(0.06) 


0.5(0.1) 


G 345.01 + 1.80 (-20,20) 


11.8(0.3) 


14.3(0.3) 


13.4(0.5) 


-14.07(0.01) 


3.08(0.02) 


0.62(0.03) 


G 351.41+0.64 (20,100) 


16.3(0.3) 


18.2(0.2) 


20.6(0.5) 


-4.20(0.01) 


3.65(0.02) 


0.86(0.03) 
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the parameters of the whole core are also given in Table|5] For 
three of them (G 265.14, G 269. 11 and G 285.26) physical pa- 
rameters have been calculated. 

3.4. Physical parameters of clumps 

In column 7 of Table|5]mean Une widths averaged over the half 
maximum intensity regions are given. For overlapping clumps 
that are not separated at this level, hne widths have been aver- 



aged over the whole regions. The following formula was used 
for weighted averages: 

(AV) . ^ , (3) 

where AV,- - velocity dispersion in the i-th position, w, - cr~'^ 
and cr, - uncertainty of velocity dispersion in the i-th position 
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where 



Fig. 2. The N2H^(l-0) spectrum toward the intensity peak of 
G 351.41+0.64. The positions of hyperfine components are in- 
dicated. The relative amplitudes of the components correspond 
to the LTE optically thin case 
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Fig. 3. Ratio of (12-12) to (23-12) peak temperatures versus 
ratio of (01-12) to (23-12) peak temperatures (Rn and Ru2, 
respectively). The LTE curve is also shown 



calculated from the fit. Uncertainties of mean line width were 
calculated as: 



o-(AV) - cr- 



a - 



(5) 



is the uncertainty of individual line width measurements taking 
into account point-to-point differences in line widths and using 
the same weights. 

Mean line widths vary from 1 km s"' to 4.7 km s"', i.e. 
much higher than thermal widths. So these values should be 
considered as non-thermal velocity dispersion in the sources 
without any recalculations. 

Following Paper I we have calculated virial masses of ho- 
mogeneous spherically-symmetrical clumps with no external 
pressure and no magnetic field (e.g. Pirogov & Zinchenko 

mm -. 



M™-(M0) = 105<Ay)2 ■ d 



(6) 



where AV and d are in km s ' and pc units. Virial masses (col- 
umn 8 of Table |5} lie in the range: ~ 30 - 3000 Mq. These 
values should be treated as upper limits. If clumps have inner 
density profiles ~ r equation Q should be multiplied by the 
factor I • which for o- > is < 1 . Mean volume densities 



3 5^2a 
5 ■ 3 

for spherically-symmetric clumps have been calculated as 



= 3Mvi,/(47rr^7i) 



(7) 



where r = d/2 is the radius of the half maximum intensity 
region, m - 2.33 amu, mean molecular mass. Mean volume 
densities lie in the range: (0.5 - 17.1) ■ lO'* cm"-^ (see column 
9 of Table|5}. The N2H^ abundances are given in column 10 of 
TablelU 



X(N2H+) 



(jV(N2H+)) 
<A^(H2)) 



(8) 



1.2 



where (A^(N2H^)) is mean column density averaged over the 
emission region defined by ///max ^ 0.5 and (A^(H2)) is 
molecular hydrogen column density calculated for the same 
region using virial mass estimate. Calculated abundances he 
in the range: (1.2 - 12.8) ■ 10"'". Averaged over 36 analyzed 
clumps, (X(N2H+)) = (5.2 + 0.5) ■ lO""'. 

4. Distributions of the observed parameters over 
the sources 

According to the present theories of star formation radial de- 
pendences of density and velocity dispersion (which is consid- 
ered to be related to magnetic field distribution) are the most 
important parameters describing the evolution of star forming 
cores. In particular, the standard theory of low-mass star forma- 
tion from a singular isothermal sphere (Shu ll977t gives a radial 
density law oc r^^ for nearly static envelope while in the logat- 
ropic sphere model for massive star formation density varies as 
r ' (McLaughlin & PudritzQ996). Velocity dispersion in the 
latter model is nearly constant in the center and rises with dis- 



tance in the outer regions as 
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Recent studies of dust continuum emission in star forming 
low-mass and massive cores have not revealed significant dif- 
ferences between density power-law indices (a) for these two 
classes of objects, implying a ~ 1.5-2. Molecular multi-line 
modeling of massive cores shows somewhat lower values (see 
Section [Q. An analysis of the N2H^(l-0) data in low-mass 
cores (Paper I) revealed the following ranges: a = 1.8 - 2.8 
for cores with stars, and a = 1 .6 - 1 .9 for cores without stars. 

Goodman et al. <1998> found that velocity dispersion pro- 
files in low-mass cores are nearly constant within so-called "co- 
herent regions" with sizes of ~ 0.1 pc, while rising in the outer 
regions with power-law index ^ 0.2 which is more in line with 
the logatropic model. Low-mass cores observed in N2H^(l-0) 
(Paper I) imply "coherent lengths" of ~ 0.01 pc but show no 
common trends in the outer regions. In contrast, molecular line 
widths in massive cores increase toward the center (Zinchenko 
[T9951 Lapinov et al. 1998 Fontani et al. 2002 1 which could be 
connected either with optical depth effects or with an enhance- 
ment of dynamical activity of gas in central regions. 

In this section we make an attempt to derive N2H^ in- 
tegrated intensity and line width radial profiles on the cores 
which, under certain assumptions, could be connected with 
density and velocity dispersion profiles, respectively. An anal- 
ysis of velocity gradients in the cores which could give a clue 
to their rotational properties is also given. 
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Fig. 4. Four examples of integrated intensity profiles together 
with the fitted power-law functions convolved with a gaussian 
telescope beam. The data lower than the dashed lines (20% of 
peak integrated intensity) were neglected 



4.1. Integrated intensity profiles 

As N2H^ emission is probably optically thin toward most of 
the map positions, except peak positions in few of them (see 
Section]^), their integrated intensity distributions should fol- 
low column density distributions if excitation conditions do not 
vary inside the sources. If in addition N2H+ abundance is con- 
stant, the integrated intensity distribution is directly related to 
the hydrogen column density distribution. Then, for infinite, 
spherically-symmetric core the power-law index for volume 
density profile will be p + 1, where p is an integrated inten- 
sity power-law index (see e.g. Motte & Andre '2001'). 

We performed fitting of N2H+(l-0) integrated intensity 
maps for 26 clumps which are not highly elongated (axial ratios 
< 2) using power-law radial distributions (oc r^i') convolved 
in 2D with a gaussian telescope beam. Although power laws 
are not optimal for fitting the radial distributions (e.g. gaussian 
functions usually give lower values) we used them in order 
to probe diff'erent star formation theories. The intensity maps 
integrated over total line range as well as those integrated over 
the (01-12) component have been analyzed. The latter should 
be less affected by optical depth effects than total integrated 
intensity maps, however, due to lower intensities these maps 
are more noisy and in several cases fits were not successful. 
Power-law index, amplitude and coordinates of maximum po- 
sition have been varied in order to get the best fit. Power-law 
indices calculated for these two kinds of maps {p and p', re- 
spectively) with fitting uncertainties are given in Table |6l We 
have also calculated power-law indices for reduced integrated 
intensity maps which contain only points with intensities equal 
or higher than 20% and 40% of peak integrated intensity. The 



calculated power-law indices (p2Q and p4o, respectively) and 
corresponding uncertainties are also given in Tabled In Fig.lH 
four examples of integrated intensity maps as ID profiles are 
shown together with fitting results for reduced maps. The ra- 
dial distance (r) in Fig.l^is the distance from the maximum of 
fitting curve to the given point of the map. The data below the 
dashed lines (20% of peak integrated intensity) have not been 
taken into account during fitting. 

A comparison of the different power-law indices derived 
for each individual source (Table |6j shows that p and p' are 
nearly the same confirming the idea that optical depth effects 
are not important in general. Mean power-law indices aver- 
aged over the sources in Table |6l are: {p} - 1.29(0.05) and 
{p'} - 1.34(0.06). We have fitted total integrated intensity 
maps without data points with intensities larger than 90% of 
the peak, which could mainly suffer from optical depth ef- 
fects. In this case the power-law indices are also close to p. 
However, a further comparison of power-law indices calculated 
for total and reduced maps reveals the following inequality: 
p > P2Q ~ P40 for most of the sources, indicating that N2H^(1- 
0) intensity distributions are steeper away from map centers 
comparing with inner regions and do not follow single power 
law. Note, that the p2Q and p4o indices calculated over the re- 
gions with relatively high intensities could be more influenced 
by possible saturation effects than the p indices. The effect of 
steepening of dust intensity distributions in outer regions was 
found by Beuther et al. (2002) for most of massive star forming 
cores they studied; they associated this phenomenon with finite 
sizes of the cores. N2H^(l-0) integrated intensity profiles in 
low mass cores also demonstrated the same effect (Paper I) but 
it was explained by fast drop in excitation temperature caused 
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Table 6. Power-law indices calculated from fits of power-law 
distributions convolved with the telescope beam to N2H^(l-0) 
integrated intensity maps with axial ratios < 2 



Source p p' p20 Pio 



RNO IB 


0.91(0.05) 


0.89(0.09) 


0.83(0.05) 


0.68(0.09) 


NGC 2B1 


1.42(0.06) 


1.34(0.16) 


1.17(0.08) 


1.19(0.02) 


S 187 (1) 






79fO 051 


n 94f0 091 


S 187 (3) 


1. 19(0.22) 




0.82(0.14) 


0.76(0.09) 


V.J I J D ^ I . jLZ^ \ l f 


q 1 f n 20^ 




1 37fO 101 




S 199 






1A(0 1 1 1 


3QC0 9Sl 


S 201 


1 3nrn i 




1 09(^0 901 


ftKCi 90^ 
u.u / \ \j.z.y ) 


Ci 1 1^ 4S-(-9 4S 

V.J 1 / J ^*-tO^ jL.*-t^ 


L .J L \ J 




1 OSfO Oft\ 




G 173.58+2.44 


1.07(0.08) 


Q7C0 


0.84(0.07) 


0.70(0. 10) 


AFPiT (=i36fi nt 




1 9nrn 

1 . L\j\\J ) 




1 97('0 961 




1.43(0.1 1) 


1 . 1 D\\J ,U\J ) 




1 91 (0 371 


S 76E 


1 99('0 03^ 


1 1 3(^0 04^ 




Q3(0 071 
\j.yj\\j.\j 1 ) 


DR21 (1) 


1.42(0.03) 


1.47(0.03) 


1.12(0.09) 


1.04(0.09) 


DR21 (2) 


1.15(0.08) 


1.14(0.09) 


0.78(0.12) 


0.51(0.20) 


S 140(1) 


0.99(0.05) 


1.02(0.08) 


0.81(0.05) 


0.65(0.08) 


S 140(2) 


1.19(0.06) 


1.11(0.10) 


0.80(0.13) 


0.54(0.14) 


S 153 (1) 


1.29(0.14) 


1.30(0.32) 


1.12(0.09) 




S 153 (2) 


1.08(0.11) 


1.23(0.25) 


0.66(0.08) 


0.46(0.08) 


G 268.42-0.85 


1.35(0.10) 


1.58(0.26) 


0.97(0.13) 


0.72(0.31) 


G 291.27-0.71 (1) 


1.33(0.06) 


1.41(0.10) 


1.05(0.05) 


0.94(0.08) 


G 345.01 + 1.80(1) 


1.35(0.03) 




0.99(0.10) 


0.69(0.15) 


G 345.01 + 1. 80 (2) 


1.37(0.03) 


1.40(0.04) 


1.05(0.06) 


0.76(0.09) 


G 345.41-0.94 (2) 


0.75(0.09) 


0.89(0.17) 


0.72(0.08) 


0.63(0.05) 


G 351.41+0.64(1) 


1.40(0.08) 


1.41(0.03) 


1.01(0.01) 


0.77(0.04) 


G 351.41+0.64 (2) 


1.30(0.02) 


1.36(0.05) 


1.01(0.03) 


0.90(0.04) 


G 351.41+0.64 (3) 


1.24(0.05) 


1.34(0.08) 


0.90(0.08) 


0.77(0.11) 



by density drop in outer core regions. The latter seems doubtful 
since the same cores mapped in CS are usually more extended 
than in N2H+ (Zinchenko et al. . 1994 1995 1998 : Juvelal996 1, 
instead, the fast decrease of intensities could be caused by 
N2H^ abundance drop (A. Lapinov, private communication). 

Mean power-law indices for reduced maps averaged over 
the sources in Table |6l are: {pia) - 1.05(0.06) and {pAo) - 
0.96(0.10). Assuming excitation conditions are constant and 
saturation effects are small in the interior regions where 
I/Imax ^ 0.2, the derived power-law indices more likely cor- 
respond to density profiles oc r"^ as the theory of isothermal 
sphere predicts (Shu 1977 1. This is in general consistent with 
mean power-law index 1.6(0.5) found by Beuther et al. (2002 1 
for radial density profiles in inner clump regions (< 32") of 
massive star forming cores as well as with the value 1.8(0.4) 
derived by Mueller et al. d2002t for dense cores associated 
with water masers. Our sample includes seven sources ex- 
plored by Mueller et al. ( 2002 1, and for five of them (RNO IB, 
NGC 281, G 173.48, S 76E and S 140) we have got intensity 
profile power-law indices. However, density power-law indices 
derived by Mueller et al. (2002) are systematically lower by 
~ 0.4 - 0.6 than the values p2o + 1 expected for density pro- 
files from our data for inner clump regions. Model calculations 
of N2H^ excitation in dense cores as well as maps of optically 
thin dust emission for the rest of our core sample are highly 
desirable to explain the origin of the discrepancy. 

4.2. Velocity dispersion profiles 

The values of line widths obtained from spectra fits are much 
larger than thermal widths (kinetic, dust or peak CO(l-O) tem- 
peratures for the sample sources are ~ 20 - 50 K, see Schreyer 
et al. [19961 Zinchenko et al. [T9951 [IWTl [T^D and could be 
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used as a measure for dispersion of non-thermal velocities in 
the sources. This parameter is widely used in theoretical mod- 
els describing non-thermal motions of gas within dense cores 
and is important for calculations of equation of state, stability 
and evolution of dense envelopes around newborn stars. 

We have compared line width and integrated intensity maps 
and found that in most cases the shapes of the regions with 
constant line width do not follow directly integrated intensity 
distributions demonstrating more complex clumpy structures. 

Since one-dimensional line width profiles demonstrate 
large scatter, we used the method from Paper 1 to calculate a de- 
pendence of averaged line widths on impact parameters. Only 
the values higher than 3(t were taken for averaging. Here, the 
impact parameter (b) is the square root of A/n, where A is the 
area of the region enclosing all the points with intensity larger 
than certain level, which in turn varies from /^ax to 0.05 ■ /^ax 
with a step 0.01 ■ /max- If the difference bt - bj-i becomes higher 
than half of the mapping step (10") the averaged line width and 
its uncertainty have been calculated for this region according to 
the formulas (|3H3. 

In Fig. |5l velocity dispersion profiles are given for ten 
sources having at least 5 points per profile including seven 
cores with one clump and three clumps separated from their 
neighbours at half maximum intensity level. The number of av- 
eraged values per point on the plots increases with b from 1 to 
10-15 at the edge. Uncertainties of averaged line widths also 
increase with b (if they are lower than (2 cr7^) where cr, 
is the line width uncertainty in the i-th map position calculated 
from the spectral fit, the latter values are shown). The cores 
with one clump demonstrate that velocity dispersion decreases 
with impact parameter or remains nearly constant in one case 
(NGC 281). In 5 of these 7 cases, IRAS sources (taking into 
account their position uncertainties) are located within 90% of 
peak intensity level. In the three remaining clumps in Fig. |5] 
there is no clear AV - b trend and either no associated IRAS 
source or it is located outside half maximum intensity level. 
RNO IB, S 199 and G 268.42 demonstrate the steepest slopes. 
Differences between central line widths and those at the edges 
of these cores are ~ 1-1.5 km s"'. For eight cores with more 
than one clump per core and at least 5 points per profile which 
are not shown in Fig. [S] there is a decrease of line width with 
impact parameter at least in two cases, G 269.1 1 and G 264.28. 
In other cases no clear trend has been found. 

An enhancement of line width in central regions in several 
cases could be at least in part due to optical depth effects. In 
fact, line optical depths ~ 0.5 - 1.5 can result in ~ 1.1-1.3 
gaussian line broadening due to saturation (e.g. Phillips et al. 
19791. Therefore, the line width profiles in RNO IB, S 76E 
(optical depths are given in Table |4} and, probably, in S 199, 
G 173.48 and G 173.58, where optical depths are more uncer- 
tain (1.3±0.7, 0.2±0.1 and 0.6±0.3, respectively), could par- 
tially suffer from this effect but cannot alone explain the ob- 
served trends in all the cases. The conclusion will not change 
if the objects consist of small unresolved clumps. 

A plausible explanation could be found in a larger degree 
of dynamical activity of gas in central regions in the vicinity of 
IRAS point sources, including differential rotation, infall mo- 
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Table 7. Results of gradient fitting 
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Fig. 5. Radial velocity dispersion profiles for the cores with one 
clump and clumps separated at half maximum intensity level 



tions and turbulence due to winds and outflows from massive 
stars. 



4.3. Velocity gradients 

It is known that an analysis of spatial distributions of line-of- 
sight velocities is fruitful in revealing rotational motions of gas 
(Goodman et al.'1993'l. In particular, solid body rotation should 
produce linear gradients on Vlsr maps. Velocity gradients also 
have been interpreted in terms of turbulent motions by Burkert 
& Bodenheimer ('5UUU'). 

In order to investigate velocity structure of the sample cores 
we have used a least squares fitting program (Goodman et 



al ll993l Paper I) which calculates local and total linear velocity 
gradients for a given core (see Paper I for a detailed explanation 
of the program). The results of fits are summarized in Table0 
where the magnitude (F) and the direction (0^, direction of 
increasing velocity, measured east of north) of the total veloc- 
ity gradient in each core, the product of the velocity gradient 
and core size (F x d), and f5, the ratio of rotational to gravi- 
tational energy (Goodman et al. 1993 Paper I) are given. We 
have not separated all the cores into clumps (as has been done 
in Section im except those having clumps that do not overlap 
at half maximum intensity level (S 187, G 133.69, AFGL 6366, 
S 255 and G 345.41). 

In Fig.|S]local velocity gradient maps (white arrows) over- 
layed over integrated intensity maps (grey scaled) are shown 
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for 20 cores. Two objects (S 187 and G 133.69) are not shown 
due to small number of points per local gradient map. The ar- 
rows point towards increasing Vlsr- Total velocity gradients 
calculated for the whole cores or distinct clumps are shown by 
bold arrows; their values are also indicated. Many sources show 
systematic and at least in some parts of maps nearly constant 
(both in magnitude and direction) velocity gradient fields, im- 
plying nearly uniform rotation (e.g. RNO IB, S 201, G 173.48, 
G 173.58, S 255, G 264.28, G 265.14). 

The cores with (not separated) closely located clumps usu- 
ally have more complex velocity fields than those composed 
of single clump or more clumps separated at the half maxi- 
mum intensity level. No direct link between velocity gradient 
distribution and IRAS source location has been found. The di- 
rection of constant velocity gradient field in several cases is 
close to elongation angle of the core. For 17 cores resolved in 
both directions there is a correlation between direction angle 
of total velocity gradient (0r) and elongation angle found from 
2D gaussian fitting of integrated intensity maps, with correla- 
tion coefficient, cc = 0.9. This fact satisfies an assumption that 
elongation of cores and clumps could be due to rotation. Yet, /3 
values vary from 4 ■ lO"'* to 7.1 ■ 10"^ with the mean value of 
0.01 for 27 sources (Table0 last column). These ratios are the 
lowest ones among dense cores with 0.2-2 pc sizes (Phillips 
.1999 ). Thus, rotation should not play significant role in core 
dynamics. 

5. Correlations between mean velocity 
dispersions, column densities and sizes 

We have searched for possible correlations between mean line 
widths (or mean non-thermal velocity dispersions) and emis- 
sion region sizes for the cores and clumps found by deconvo- 
lution method (see Section l33i . Beginning from the works of 
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KCC2ai 




Fig. 6. Local velocity gradient maps superposed on integrated 
intensity maps. Total velocity gradients are shown by bold ar- 
rows together with their values (in km s ' pc '). IRAS point 
sources are marked by stars. The uncertainty ellipses are also 
shown 




Fig. 6. Continued 



Larson ([TWT) and Myers (11985> such correlations are widely 
used in estimates of clouds stability and turbulence. 

In Fig. 0^a) we plot mean line widths versus emission re- 
gion sizes in logarithmic scale for 26 cores and clumps having 
size estimates (Table |5|i, with the exception of three cores pre- 
senting more than one velocity component (G 3 16.77, G 345.01 
and G 351.41). The correlation coefficient between parame- 




Fig. 6. Continued 



ters is not high (~ 0.5), yet, the correlation is significant. The 
hypothesis of zero correlation is rejected at 0.6% significance 
level (Press et al. ll992> . The slope of the regression line, cal- 
culated by the ordinary least-squares method, is 0.3(0.1), and 
the intercept coefficient is 0.38(0.03). Practically, the same re- 
sults have been obtained excluding the cores with steepest line 
width profiles (RNO IB, S 199 and G 268.42, see SectionlO. 
for which mean line widths could be enhanced due to additional 
broadening in the central regions of the maps. However, since 
the compared variables are independent, it could be more cor- 
rect to use linear regression methods that treat them symmetri- 
cally (e.g. orthogonal, reduced major-axis or bisector methods, 
see Isobe et al. l990). Applying these methods to the data we 
have obtained higher slopes (0.4-0.65) and nearly the same in- 
tercept values. 

The regression line calculated by the ordinary least-squares 
method is shown in Fig. 0a) as a solid line together with two 
additional linear regressions found by other authors with the 
same method for massive cores. Whereas the dependence of 
Caselli & Myers ( 1 995) (dashed line) was calculated using 
multi-line data for more compact cores in Orion A and B, the 
dependence of Zinchenko f2000'^ (dotted line) was calculated 
for the sample observed in CS and C^'*S. All the slopes agree 
within the uncertainties of their determination. 

We have also searched for possible correlations between 
peak N2H^ column densities (which in the case of con- 
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Fig. 6. Continued 

stant N2H^ abundances trace hydrogen column densities) and 
sizes (Fig. Elb)). These parameters have been chosen be- 
cause they are independent, unlike to mean densities and virial 
masses which have been calculated using size estimates (see 
Section]^}. No significant coiTelation has been found. No cor- 
relation has been also found between mean N2H^ column den- 
sities (averaged over emission regions with ///max ^ 0.5) and 
sizes. 

It is doubtful that the observed (AV) - d dependence is con- 
nected with general turbulence properties or with some com- 
mon equation of state for the sample cores. As it is shown 
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Fig. 7. Mean line widths versus emission region sizes (a) to- 
gether with regression (solid) line. Dashed and dotted lines cor- 
respond to the line width - size relations from Caselli & Myers 
(Q,995 ) and Zinchenko ( 2000), respectively. Peak N2H^ column 



3 densities versus emission region sizes (b) 



in Section 14.21 there is no evidence for growing line width 
with distance inside the cores. On the other hand, the de- 
pendence found does not contradict the assumption of grav- 
itationally bounded objects in virial equilibrium. Depending 
on the so-called virial parameter, MyiJM (Bertoldi & McKee 
1 19921 1. where virial mass is calculated from equation (|6jl and 
M is an independent estimate of mass, such objects could be 
in the critical state on the threshold of gravitational collapse 
(MyiJM ~l) or in pressure equilibrium (Mvk/M >>1) (e.g. 
Pirogov & Zinchenko'1998"(. If the cores are close to the crit- 
ical state, with magnetic energy being of the order of kinetic 
energy, the slope of the (AV) - d dependence should be equal 
to 0.5, with intercept coefficient depending on external pressure 
and inner density profile. Column density in this case does not 
depend on size. If the cores are in pressure equilibrium, any 
(AV) cc d'' dependence implies a corresponding M oc d^^^'' 
dependence. 

Our virial masses for four cores coincide within a factor of 
1.5 with the M„ masses found by Mueller et al. f2002"), and in 
one case (RNO IB) the virial mass is ~ 5 times lower than the 
corresponding M„ value, which could be connected with dif- 
ferences in core size. Masses derived from CS observations for 
most of the sample sources (Zinchenko et al. . 1994.. 1995. .19981 
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Juvela [T996l are several times larger than our virial masses (ex- 
cept G 316.77), yet, CS emission is usually more extended. 
Therefore, if the studied cores are in equilibrium, they could 
be close to the critical state rather than in pressure equilibrium. 

6. Comparison with physical parameters of 
low-mass cores 

In Table|8lphysical parameters of the clumps averaged over the 
studied sample are given together with the same parameters 
derived in Paper I for low-mass cores with and without stars. 
Physical parameters for massive cores are taken from Table |3l 
and Table |5l We have used weighted averages for the param- 
eters which uncertainties resulted from fits (line width, size, 
axial ratio, power-law index, velocity gradient), for the others 
standard arithmetic means are given. In brackets we give stan- 
dard deviations (as in Paper I). 

Apart from well-known differences in virial mass and ve- 
locity dispersion between these two samples, one can see that 
massive cores have higher sizes and column densities, lower 
mean densities and velocity gradients than low-mass counter- 
parts. Note, that for low-mass cores d - 2r, where r was taken 
from Table 6 of Paper I. Mean power-law indices for integrated 
intensity maps for the studied cores are also given as well as the 
same parameters for low-mass cores (see Table 6 from Paper 1, 
p - a - 1). On average, integrated intensity profiles in massive 
cores have nearly the same power-law indices as in low-mass 
cores within uncertainties of their determination. Mean axial 
and p ratios, as well as mean N2H^ abundances also are nearly 
the same taking into account their uncertainties. Velocity dis- 
persion in massive cores is either constant or decreases with 
distance from the center in contrast to low-mass cores where 
no common tendency was found. 

The fact that mean density in massive cores is lower than in 
low-mass cores is somewhat unexpected. For example, Fontani 
et al. (|2002i recently found densities of 10^ - 10^ cm"^ for 
the sample of massive cores associated with ultracompact H II 
regions and having sizes 0.2-1.6 pc. Possible explanation could 
be found if dense cores from our sample actually consist of 
unresolved clumps much smaller than the telescope beam. If 
density of such unresolved clumps is the same as for low-mass 
cores with stars (see Table |8} the clumps volume filUng factor 
should be ~ 0.3. 

7. Conclusions 

In order to derive physical conditions of dense gas in molecular 
cloud cores where massive stars and star clusters are formed, 
35 objects from northern and southern hemispheres have been 
observed in the N2H^(l-0) line. N2H^ emission has been de- 
tected in 33 sources, and detailed maps have been obtained for 
28 sources. The results can be summarized as follows: 

1. PeakLTEN2H^ column densities lie in the range: 3.610'^- 
1.5- 10'** cm.-^. Intensity ratios of (01-12) to (23-12) hyper- 
fine components in most cases are slightly higher than the 
LTE value. The optical depth of (23-12) component toward 
peak intensity positions of 10 sources is ~ 0.2 - 1. 



2. N2H^ intensity maps in most cases show elongated or 
more complex structures. Two-dimensional gaussian fit- 
tings have revealed 47 clumps in 26 objects. Their sizes 
lie in the range: 0.3-2.1 pc, the range of virial masses is 
~ 30 - 3000 Mq. Mean N2H+ abundance for 36 clumps is 
(5.2 + 0.5)- 10-1°. 

3. Power-law indices for radial integrated intensity profiles in 
the cores with axial ratios < 2 have been calculated. The 
mean value of power-law index is close to 1 .3 for total inte- 
grated intensity maps in 25 objects, and about unity for re- 
duced maps, where positions of low intensity were rejected. 
This corresponds to a density profile oc in the core inner 
regions, assuming constant N2H^ excitation conditions and 
abundances, and small saturation effects. 

4. In the cores with relatively extensive and high quality maps 
N2H^ line widths either decrease or stay constant with dis- 
tance from the center implying an enhanced dynamical ac- 
tivity in the central regions in the vicinity of IRAS sources. 

5. Elongation angles of the 17 cores and distinct clumps cor- 
relate with total velocity gradient direction confirming that 
elongation could be due to rotation. However, the ratio of 
rotational to gravitational energy (4 ■ 10^^-7.1 ■ 10"^) is too 
low for rotation to play a significant role in the dynamics of 
the cores. 

6. A correlation has been found between mean line widths 
and sizes of the clumps, (AV) oc £/"^=^'^ '. No correlation 
has been found between peak N2H^ column densities and 
the sizes. These results do not contradict the assumption of 
gravitationally bounded objects in virial equilibrium. 

7. A comparison between mean physical parameters of mas- 
sive and low-mass cores (Paper I) reveals that massive cores 
have higher sizes, N2H^ column densities, virial masses 
and velocity dispersions as well as lower mean densities 
and velocity gradients than low-mass counterparts. Mean 
N2H^ abundances, axial ratios and ratios of rotational to 
gravitational energy are nearly the same. N2H^ integrated 
intensity profiles have nearly the same mean power-law in- 
dices while velocity dispersions in massive cores have a 
tendency to decrease with distance in contrast to low-mass 
cores, where no common trend was found. 
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Table 5. Physical parameters of N2H+ clumps 
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Parameter 


Massive cores 




Low-mass cores 










With stars 


Without stars 




Mean 


Number 


Mean 


Number 


Mean 


Number 


AV(kms-') 




^1 

D i 






fl Q/TJ 1 'l 


ZD 


yV(N2H+) (10'2 cm-2) 




JO 






OKi) 


LJ 


rf(pc) 






0.14(0.03) 


22 


0.1(0.02) 


13 


Axial ratio 


1.4(0.4) 


47 


1.9(0.7) 


35 


2(1) 


19 


nv,r (lO-* cm-3) 


6(4) 


36 


20(20) 


34 


12(7) 


19 


Mvir (Mo) 


658(604) 


36 


9(16) 


34 


3(2) 


19 


X(N2H+) (10-'°) 


5(3) 


36 


3(2) 


34 


2(1) 


18 




1.3(0.1) 


25 


1.2(0.1) 


9 


0.8(0.04) 


8 


r (km s"' pc"') 


0.5(0.4) 


27 


2(1) 


14 


2(1) 


12 


P 


0.01(0.02) 


27 


0.02(0.02) 


13 


0.02(0.02) 


7 



